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Abstract 



(N 
> 

We find the general solution for the spacetimes describing the interior of static black 
holes with an equation of state of the type Tq = T\ (T being the stress-energy tensor). 
This form is the one expected from taking into account different quantum effects as- 



sociated with strong gravitational fields. We recover all the particular examples found 
in the literature. We remark that all the solutions found follow the natural scheme of 
an interior core linked smoothly with the exterior solution by a transient region. We 
£jT)! also discuss their local energy properties and give the main ideas involved in a possible 

generalization of the scheme, in order to include other realistic types of sources. 

X 

1 Introduction 

Any static black hole (BH) arises from the gravitational collapse of some object. Under the 
premises in this work, the object has not (yet) shrinked indefinitely and has not given rise 
to a spacetime singularity. It is then natural to consider two regions: one exterior to the 
object, and the object itself. The exterior region, as is well known, can be described by a 
spacetime belonging to the Reissner-Nordstrom |], |2j solution. In the absence of electric 
— or magnetic — charge it is simply Schwarzschild's spacetime. Furthermore, one can add 
a cosmological constant, following recent observational results |3|, M. Then, the spacetime 
belongs to the Kottler-Trefftz solution family ||. In this case, the global properties of the 
spacetime clearly change, e.g., the spacetime is no longer asymptotically flat (see e.g. |], 01). 

Going a step further, we consider the body itself as composed of two main regions. One 
is its surface and the other the rest of the body, i.e. the interior region. One may expect 
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that some mechanism — having to do e.g. with quantum gravity — will be able to stop the 
collapse of the body. Therefore, we will think of the interior of the body as being described 
by some spacetime product of the present knowledge available on the merge of quantum field 
theory and gravitation. A widely studied issue in this direction is that of quantum vacuum 
effects HHpJl] and the resulting spacetime turns out to be a de Sitter (dS) or anti de Sitter 



(AdS) one |15|]-|23|]. There are other alternatives, as e.g., those of |25|, |26| . 

In all these cases considered so far, either no distinction has been made between the inte- 
rior and the exterior of the body (see e.g., p0|-p6|) or there appears a singular distribution 
of matter at the surface of the body (see e.g. [18|, This distribution is singular in the 
sense that it is a matter surface density — called singular shell. However, contrary to the 
case of electromagnetic charge densities, a matter surface density has neither been observed, 
nor is it predicted by any theory. It is thus more natural to assume that the matter on 
the surface of the body is distributed across the body, and leave for a subsequent study the 
issue of whether this region is thin or thick, in comparison with the region dominated by 
quantum vacuum effects, through one of the solutions referred to before. Finally, the only 
work considering all the features of the structure of a regular static BH with a clear physical 
source is |27|. However, Nariai spacetime was absent, as well as an implementation of previ- 
ous attempts and a complete study of (local) energy conditions. Thus, in our opinion, it is 
worth carrying out a unification of the different results obtained so far, as well as extending 
them in order to cover some impotant issues that were overlooked in those analysis. Here 
we provide, for the first time to our knowledge, the general solution of the scheme discussed 
above. In particular, we carry out an implementation of all those previous works which, for 
one or another reason did not comply with all the requirements already specified. We also 
perform a study of the local energy conditions in all these cases. 

Finally, it is also important to introduce other kind of solutions for the interior region, 
aside from the ones referred to before, which arise from results, or just hints, coming from 
the contribution of the quantum vacuum to gravity. To summarize, these are the points that 
will be dealt with, successively, in the body of this paper comprising the next 10 sections. 
They are clearly identified by their titles and will need no further specification here. Sect. 12 
is devoted to some final remarks, and in Sect. Owe provide the conclusions of the work. A 



brief survey can be found in 28 



Throughout this work we will use units such that G — 1, c = 1, Einstein's equations 
are written in the form G a p = 8irT a p, where G a p is the Einstein tensor — we follow the 
conventions of [jl] — and T a p is the energy- momentum tensor. A prime will denote derivation 
with respect to the coordinate r. 
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2 Spacetimes with a SSQV as a source 



Due to the imposed limitation of non-accepting singular mass shells, the spacetimes describ- 
ing the interior of the body are not allowed to be of some well known kind as dS or AdS 
spacetimes. Indeed, the interior solution cannot be everywhere a spatial isotropic solution, 
as dS or AdS, because this would lead to a sudden change in the pressures exerted by the 
body to the exterior, and would lead to the appearance of a singular mass shell f3(|. Now, 
the type of generalization depends on the underlying physics one is able to assume. The 
isotropic case is suitable in order to link it with the expected contributions of a dominating 
quantum vacuum, especially those associated with vacuum polarization. As we are dealing 
with spacetimes which are spherically symmetric, a natural generalization is to assume that 
the body may be described by a solution which is invariant to any non-rotating observer, 
with a free radial motion, instead of a solution which is invariant to any observer. This gen- 
eralization of the energy-matter content of the body is called spherically symmetric quantum 
vacuum (SSQV), after [^C§ — see also [|2[] — and requires the imposition of Tq = Tj, for 
any non-rotating observer. This is the type of enery-matter content that is considered in 



|p!8[-p9| and will be the one used in the first part of this work, until we get to Sect. [10]. In 
particular, SSQVs and non-linear electrodynamics have given some relevant results on the 
issue of regular BHs, see e.g. P3|-P5||, pE] . 

We shall now characterize the families of spacetimes that are suitable to become SSQVs. 
Any static, spherically symmetric spacetime can be conveniently described by 



ds 2 = -F(r) dt 2 + F~\r) dr 2 + G 2 (r) dtf, 



(1) 



where dfl 2 = dd 2 + sin 2 6 dip 2 . There are certainly other ways to represent these spacetimes 
— which avoid the problems occurring near the possible horizons — or by putting R 2 dQ 2 , 
provided G' = dG{r)/dr ^ (see e.g. f3(J, [31]]). For a local observer at rest with respect to 
the coordinate grid of ([[]), a standard computation of T Q/ g yields (p is the energy density, p 
the radial pressure and p 2 , P3 the tangential pressures, measured by this observer) 



1 

G 2 



- F(G 
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8np 2 = 8np 3 = 

Imposing p + p = in Eqs. @, (|3|), we get 



+ 2GG") - GG'F' 

1 + FG' 2 + GG'F' 

F" FG" F'G' 
~2~ + ~G~ + ~G~' 



FG" = 0. 



(2) 
(3) 
(4) 

(5) 



We now use that G cannot be zero in any open region. Two alternatives appear: F = or 
G" = 0. If F = 0, the expression ([[]) is useless. It is first necessary to change the coordinate 
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system of ([!]) by dT = dt + (1 — F)/Fdr, while keeping the rest unchanged. Then one can 
impose F = 0. The result is 



ds 2 = 2dTdr + 2dr 2 + G 2 {r)dQ 2 . (6) 

In the orthonormalized cobasis given by 0° = dT/\r2, 1 = dT/y/2 + y/2dr, 2 = Gd6, 
3 = G sin 6 dcp, the Ricci tensor takes the form 

Ricci = ^(-0° ® 0° + 0° ® 1 + 1 ® 0° - 1 ® 1 ) 

G l (7) 
+ ^(@ 2 ® 2 + @ 3 ®0 3 ). 

On the other hand, for a SSQV we must have p + p = (the conditions on T a p being directly 
translated into conditions for R a p) 

Ricci = R 00 (-®° N ® ®° N + 0jv ® 0jv) + i? 22 (©^ 0^ + 0^® 0|r), (8) 

where {0^} is some orthonormalized cobasis, not necessarily coincident with the one used 
in the computation of (0). Therefore, we must look for an orthonormalized cobasis for 
which the Ricci tensor (|7]) becomes of the type (|8|). Clearly this is the same as finding out 
whether we can have linear expressions ®° N = A@° + B® 1 , and ® N = C®° + D® 1 , with 
— ®° N -®° N = ©^-©jy = ©^-©jy + l = 1. However, Eq. @ is invariant under these changes. 
The only solution that makes (0) and @ compatible is then 

G" = 0. (9) 

If F ^ 0, we also have G" = 0. Thus G" = constitutes the proper characterization of any 
possibility. 

Now, from G" = two distinct alternatives appear 

G = j, or G = ar + 7, (10) 

where a(^ 0) and 7 are constant. Only the latter has been considered in detail in the 
literature of regular BHs. We will study it in the sequel. 



2.1 Other expressions for the spacetimes describing SSQVs 

In order to include the possible horizons, we write the metrics ([!]) under the common form 

ds 2 = -(l-H) dT 2 + 2HdTdr + (1 + H) dr 2 + 7 2 rffi 2 , (11) 
ds 2 = -(l-H)dT 2 + 2HdTdr + {l + H)dr 2 + r 2 dQ 2 , (12) 

where H = 1 — F, and the coordinate change is given by dT = dt + (1 — F)/ Fdr. We will use 
these forms in the sequel. We have also used the fact that the case g = ar + 7 is physically 
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equivalent to the case G = r. This is intuitively seen because a merely represents the scale 
of units used for r and 7 is an arbitrary (constant) origin. In terms of coordinate changes we 
have: The metric (j^) for G = ar+7 is ds 2 = —F(r) dt 2 -\- F~ x {r) dr 2 +(ar+ / y) 2 dQ 2 . Recalling 
that a^O, one can define a new radial coordinate f = ar+7. The metric becomes then ds 2 = 
— F[(r — 7) /a] dt 2 + a~ 2 F~ l [(r — 7)/a] <if 2 + f 2 rffi 2 . Now, under a reparametrization of the t 
coordinate by <it = adt we get (is 2 = — a 2 F[(r — j)/a] dt 2 + a~ 2 F~ l [{r — j)/a] df 2 + f 2 dfi 2 . 
Whence one can conclude that any member of (|l|) with G = or +7 with a 7^ is equivalent to 
another member of (P with a = 1 and 7 = 0. Since we are studying the general description 
of SSQVs it is enough to consider the representation a = 1, 7 = and arbitrary F(r) to 
include any case of SSQV. 



Furthermore, it is easy to show that (12) can be written in the Kerr-Schild form [32|: 

ds 2 = ds 2 + 2H(r)£ <g> £, (13) 

where ds 2 stands for the flat spacetime metric, H is an arbitrary function of r and £ is a 
geodesic radial null one-form, £ = (l/^/2)(dt ± dr). Thus, the SSQVs in (|TJ) can be thought 
as the family of maximal spherically symmetric spacetimes expanded by a geodesic radial 
null one-form from flat spacetime (GRNSS spaces). 

To summarize, there are only two — non-equivalent — families of SSQV. The case with 
G 1 = is characteristic of the Nariai solution p3| , |34|. The Nariai solution is a solution 
of Einstein's equations for the same pattern as the de Sitter solution, i.e. T a/ g = Aog a p, 
being A the cosmological constant. The difference lies in the "radial" coordinate. In the 
Nariai case there is no proper center for the spherical symmetry. Therefore, we shall call the 
spacetimes with G' = generalized Nariai metrics. Finally, the other case corresponds to 
the GNRSS which constitute a distinguished family of the class of Kerr-Schild metrics. 



3 Geometrical properties of the solutions 

3.1 Generalized Nariai metrics 

Using an orthonormal cobasis defined as @° = (1 — -y) dT — ^ dr, 1 = (l + y)<ir+y dT, 
2 = 7 d6, 3 = 7 sin 9 dip, we see that the Riemann tensor has as independent components 

#0101 = ~H"/2, #2323 = 1/7 2 - (14) 

The Ricci tensor is characterized by Rqo = —Rn = —H"/2, R22 = R33 = 1/7 2 - The scalar 
curvature is R = H" + 2/7 2 , and the Einstein tensor has the following non-zero components 

Goo = -Gu = I/7 2 , G 22 = G 33 = -H"/2, (15) 

The isotropic solution — the one to be found at the core — yields H = (l/7 2 )r 2 + br + c, 
where b and c are arbitrary constants. Without losing generality, we can set b, c = (as 
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they are clearly gauge freedoms for any spacetime in the family). Thanks to the presence of 
the Nariai solution inside this family — T^ anai = A g a p — the factor I/7 2 can be identified 
with Aq. Thus, the only isotropic quantum vacuum belonging to this family is the Nariai 
solution. 



3.2 The GNRSS metrics 

First we note that these spacetimes fulfill the relation 

P2 = -P- ^f- (16) 

As a consequence, for a regular source p 2 — > —p as r = is approached. Therefore, in 
any regular solution the spacetime becomes more and more isotropic as r — ► 0. Thus, the 
contribution from the quantum vacuum becomes more and more dominant as r tends to 
and the initial idea of distributing the singular mass shell across some region of the body is 
completed. 

On the other hand, one can choose a similar cobasis as in the Nariai-like case, just 
replacing 7 by r. The Riemann tensor has the following independent components, 

R - H " R -R - H ' 

"0101 — TTi -n-0202 — -n-0303 — — 7~; 

2 2r 

H' H 

-R1212 — -R1313 — 7; ? -R2323 — (17) 
2r r z 

The Ricci tensor for these spacetimes has the following non-zero components Rqo = —Ru — 
-{l/2)[H" + (2H'/r)], R22 = R33 = 0-/r)[H' + {H/r)\. And the scalar curvature is given 
by R = H" + (AH'/r) + (2H/r 2 ). Finally, the non-zero components of the Einstein's are 
Goo = = {l/r)[H' + (H/r)), G 22 = G 33 = -(1/2) + (2H'/r)}. Other expressions 

that will be used later are 

1 G' r 
Coo = —G\i = —(Hr)', G 2 2 = C 33 = G o 5 — • (18) 

In this case, the isotropic (regular) GNRSS is the de Sitter solution, given by H{r) = 
(A /3)r 2 . 



3.2.1 Exterior metrics and GNRSS metrics 

It turns out that all of the possible exterior metrics — see Sect. |I] — also belong to the 
GNRSS family. The function H is H(r) = (A ext /3)r 2 + 2m/ r — q 2 /r 2 , where A ext stands for 
the external cosmological constant, m is the ADM mass of the BH and q its electromagnetic 
charge. 

This coincidence will be very useful in the following section. 



6 



4 Junction of the interior and exterior solutions 



The junction, or matching, of two spherically symmetric spacetimes is well-known (see e.g. 
||36|-|]tjy|). The general form of a hypersurface that clearly adjusts itself to the spherical 
symmetry of any of these spacetimes is as follows 

• e = a,, 

E: r m (19) 

r = r[X), 
J = t(X), 

where {X,Xe,X v } are the parameters of the hypersurface. One must thus identify both 
hypersurfaces in some way. The identification of (Aj)i with (X{) 2 (1 and 2 label each of the 
spacetimes) is the most natural one, due to the symmetry of the above scheme. In the sequel, 
1 labels the exterior spacetime and 2 the interior one. 

In order to match the exterior solution with the interior one, one basically demands the 
coincidence of the first and second fundamental forms of S at each spacetime — the other 
way is to accept the presence of singular mass shells, which would not require the coincidence 
of the second fundamental form, but we will dismiss such unphysical option. In order to 
include the possibility of matching the interior and the exterior at null hypersurfaces (e.g. 
at an event horizon) one can follow the formalism in |[39|| . It is worth recalling that the 
exterior region is described by a member of the GNRSS family. Thus, we have to consider 
two possibilities: matching a generalized Nariai metric with a GNRSS one and two GNRSS 
metrics with each other. 

4.1 The junction of a generalized Nariai metric and a GNRSS one 



In this case one easily gets r x (A) = 7 = const., £i(A) = const. — see for full details. 
This result tells us that the junction between a generalized Nariai spacetime and a member 
of the GNRSS family is impossible. It would only happen for a (two-dimensional) surface. 
Therefore, any member of the Nariai class cannot be regarded as a good candidate in order 
to represent the interior structure of a regular, static BH. 

4.2 The junction of two GNRSS spacetimes 

In this case one gets that two members of the GNRSS family match with each other if 



and only if — see e.g. — either ri(A) = r 2 (A) = R = const., t\ = i 2 , [H] = [H'\ = 
or 7*1 + t\ = r 2 + t 2 = const. The last condition, however, describes the motion of a null 
hypersurface and is not an acceptable solution in order to describe the matter inside a static 
BH. Therefore, we reach the conclusion that: The only acceptable hypersurfaces fulfilling the 
matching conditions, that preserve the spherical symmetry, between two spacetimes of the 
GNRSS family, are those satisfying ri(A) = r 2 (A) = R = const., ii = i 2 , [H] = \H'} = 0. 
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Without losing generality, one can choose t\ — t 2 — A, because of the global existence of 
the Killing vector d t . Moreover, we realize that the chosen coordinates are privileged ones, in 
which the matching is explicitly C l . The hypersurface, E, will be timelike, null or spacelike 
according toiJ<l,iJ = lorif>l, respectively. 

To summarize, if vacuum polarization is to be the dominant quantum effect, the most 
simple way to construct a regular BH is to build it upon GNRSS spacetimes. 



5 Regular interiors of the GNRSS type 

As mentioned elsewhere, the exterior region can appropriately be characterized by a member 
of the Kotller-Trefftz class, which is a subclass of the GNRSS family. Then, the matching 
conditions between the exterior and the interior regions are: 

H 2 {R) = H 1 {R) = *£-^ + ±ie, (20) 
H' 2 {R) = H[{R) = + % + 2 -^R. (21) 

Moreover, the aim here is to focus on those interior solutions which are everywhere regular. 
From the expressions of the Riemann tensor and the metric, one sees that this may only be 
accomplished if 

H 2 (0) = 0, #£(0) = 0. (22) 

Thus, we finally encounter four conditions in order to have a regular interior solution. 

From now on, we will consider H 2 to be a smooth function of the variable f = rj R, a 
most natural hypothesis in view of the regular character prescribed for the interior solution. 
In this case, the origin conditions tell us that 

oo 

H 2 (f) = Y,b n f n . (23) 

n=2 

Now, one has to impose the two other conditions. Obviously it is the same to consider H 2 (r) 
or H 2 {r — 1) in the whole procedure. However, we will first work with H 2 {r — 1) in order to 
implement the junction conditions directly. From the preceding result, one immediately has 

oo 

H 2 = J2a n (?-lT, (24) 



n=0 



and the junction conditions tell us that 



2m q 2 : Ai , 



ao = Hi{1) = ii-h + Y R > (25) 
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where Hi(r) = (2m/R)(l/r) — (q 2 /R 2 )(l/f 2 ) + (AiR 2 /3)r 2 and a dot denotes derivation with 
respect to f. The following step is to impose regularity of the solution, Eqs. (0). We get 

oo oo 

y £(-l) n a n = 0, E(-l) n na n = 0, (27) 

n=0 n=l 

which, by virtue of the matching conditions, yield 

±(-1^= + % + (28) 

n=0 R R 6 

E(-l)> + 2 K + 2 = 2(-^ + ^ + . (29) 

It is clear that there are infinitely many possible candidates for these interiors. 

5.1 Isotropization 

Let us further analyze how they behave near the origin. Taking into account the expression 
of if 2 in powers of f and using Eqs. (|18D , we get 



Gn = -^E(^ + 1)^" 2 > G2 2 = ^ 2 t( l+ 2 %- 1 - 2 . (30) 

/ — 2 t — 2 \ / 

It is then clear that Gn and G22 are different from each other. Yet we have the very relevant 
property that, for any of these spacetimes, it holds 

limGn = lim(-Goo) = limG 22 = hmG 33 = — —. (31) 

r— i-O r^O r— >0 r— >0 /£ z 

Whence, we see that a general isotropization of the Einstein tensor — and consequently of 
the energy-momentum one — independent of the model is actually accomplished. In terms of 
a\ we get 

1 00 

G n = "I* E A Mf M , (32) 



R 2 

11 M=0 



oo / 1 

A M = (-l) M (M + 3) £ (-W ' 



v / — 2 — M ; 



and 



So that 



G 22 = ~4 E ^4,f tf (33) 

M=0 Z 



limGn = lim(-Goo) = limG 2 2 = hmG 33 = -— |. 



4, = 3E(-i) , (^ 2 Jfli. 
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Finally, making no further assumptions on the coefficients of H2, we can isolate two of 
them in terms of the rest. For simplicity, we shall isolate a 2 and 03. The result is 



a 2 



R 2 



1=4 



With this in hand we can write, respectively, the previous expression for the central value 
in terms of a\ or bi, I > 4, 



Gn(0) = G 2a (0) 



3 


8m 






3 


8m 


R 2 





+ ^ 2 + E(^-3)^ 

d Z=4 



6 Examples 

We will consider six examples. Two constitute the well-known proposals of |Tj|, [H], |4[] and 
pO], |2T|| . Two more come from a proposal in pB|-||26||, for electrically charged bodies, and 
the proposal given in p9 |, for magnetically charged ones. More specifically, we will here 
derive their corresponding analogues, within the present scheme (what is actually more than 
simply re- writing those cases). The remaining two examples constitute a family of brand 
new candidates, which naturally arise from the preceding expressions. We will start with 
this last pair. 



6.1 Two arbitrary powers 



Let us just make the choice that only two specific powers of H(r), say M and N, be present. 
In order to fulfill the regularity conditions, both must satisfy M, N >2. However, if we wish 
to obtain a de Sitter-like behavior at, and near, the origin, we must necessarily impose that 
one — and only one — of them, say M, be equal to 2. Thus, H2(f) reads -£^2(7) = ^f 2 + bNr N , 
for N > 2, with, 



N+l \ 
N-2J 

Gn(f) and G 22(f) read 




Gu(r) 



-Ai 



q 2 ( N + 2 \ 
R 2 \N-2J 

-recall Eqs. (|T 

6m (N + 1\ 



A t R 2 



(N-2)R 



—3m + 



2f 



-1) + 



i? 4 



^ 1 ~\ a 6m f N 
G 22 (r) = -A x + — 




R 4 



3{N + 2) - A(N + l)r N + 2 ' 
iV-2 

3(N + 2)-2N(N + l)r N+2 
N-2 
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Whence, one readily sees that their finite value at the origin coincides, as expected, 

r, ,nN ,nN a 6m/iV+l\ 3q 2 f N + 2\ 

Gu(0) = G 22 (0) = -A x - — — — + 4- — — . (34) 



R 3 \N -2 R 4 \N-2 



6.1.1 Lowest powers 



This example corresponds to the case in which H 2 is a polynomial of lowest degree. This 
amounts to setting a\ = 0, I > 4, in the general expressions. Its interest lies in its being the 
simplest possible situation. The result is 

rr /-n f 8m W ^iR 2 \~2 ( 6m Aq 2 \ 3 

and 

24m 15g 2 8 / 2g 2 \ 

= RF ~R^ ~R? \ m ~R~ ) r ' (36) 

G 22 = -^ + l^-A 1 + ^L-^V~ (37) 

Notice that Gn tends to — Ai — q 2 /R 4 as f tends to 1, the same value as G^'(f = 1), in 
accordance with Israel's conditions 1061. 



6.2 Israel and Poisson's model 

In reference [18] — see also [|19) — a plausible candidate for the energy- matter content of the 



interiors of regular non-charged BHs was proposed. The authors proposed that a singular 
layer of non-inflationary material should exist between the de Sitter core and the external 
Schwarzschild metric. However the usual spirit of matching a stellar interior with a vacuum 
exterior was lost, the reason being the unavoidable presence of a singular layer acting as a 



matter surface density. Indeed, in [^4j it was argued that their approach could be improved 
by imposing a smooth transition from the hypersurface to the de Sitter core. Yet this step 
was not implemented. In any case, it was the only available candidate to continue the studies 
of quantum regular BHs at that time. The task here will be to see whether this geometrical 
and physical model can be recovered from our analysis. 

In order to to do that, we search for a solution within our family which be as close 
as possible to this particular solution. What amounts to looking for a de Sitter core for 
small values of f and a quantum contribution of the type of the square of the characteristic 
curvature of Schwarzschild spacetime near the matching hypersurface. These features taken 
into account, we set for the interior 

(Goo)int = -(Gll)int = ^ - Bj .3j2> ( 38 ) 
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where A and B are two constants, to be determined. Using Eqs. (|18D, we obtain 



,2 



H 2 (r) = -r-—^i (39) 



3A(A + Br 3 ) 

where we have imposed Eqs. fl2~2"|). The matching conditions lead to 



' ^(R), l A A ■ l R " = Ih[(R), (40) 



A{A + Br 3 ) R 2 iy n A(A + BR 3 ) 2 R 

where H\ comes, as usual, from the external model. 

In the exterior region, close to the matching hypersurface, the quantum contributions do 
not turn into a cosmological-like term. They are of the form (Goo)ext. oc m 2 /r 6 , as mentioned 
before. We thus have a quantum exterior which is different from the one encountered in the 
rest of the examples and sections before, which cannot be described by a member of the 
Kottler-Trefftz class. Fortunately, our preceding results are still useful. In fact, one realizes 
that it is possible to select a suitable exterior with a similar form as (^), just by setting 
^4ext. = 0, and B~^ t = am, where a = f3L P [, being (3 of order unity, and L P1 the Planck 
length (a is of order unity in Planckian units). j3 2 is related with the number and type of 
the quantized fields, 0, This choice yields 

2m 1 fam\ 2 
Hl{r) = ~ ~ 3\—) • (41) 
where we have taken into account that the exterior region is dominated by the Schwarzschild 
geometry — with mass m — for large values of f. 
Now, using Eqs. (|40]) and ((41]), A and B yield 

/ \ 
1 



A=^-^, B 



am am 

o 



a 2 m 



(42) 



R 3 V R 3 J 

Finally, using Eqs. (j38l), (G22)int = (2Br 3 — A) /(A + Br 3 ) 2 , where A, B have been given 
above. At the origin 

G u (0) = G 22 (0) = ~. (43) 

To summarize, we have proven here that a spacetime model within our family satisfies 
all the required geometrical assumptions, and yields the particular form of Goo ; both for 
the interior and the exterior of the body, as in the above mentioned references. A more 



throughout comparison of that model and ours will be given in Sect. 1772. 



6.3 Dymnikova's model 

Some time after the appearance of the previous cases a new model for a regular interior of a 
non-charged BH was proposed |2(J. However the approach was now quite different to that of 
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the previous authors. Now Schwarzschild's solution was only recovered in an asymptotical 
sense, for f approaching infinity only. However, if a sufficiently quick convergent matter 
model was obtained, then the quantity of mass outside the horizon of the collapsed body 
could become as negligible as desired with regard to the interior mass. Thus one would, 
at least, recover a trial model, interesting enough to support or reject the conclusions of 
the previous authors. In a later work [[H]], the model was extended to incorporate the 
observational fact in favor of a non-vanishing cosmological term in the exterior region. We 
will deal in this subsection with such model, but considering a definite end to the collapsed 
body. 

The imposition for the energy-matter content for the interior will be of the form 

(Goo)int = -(Client = ^exp (-f 3 ) + B, (44) 

where A and B are two constants to be determined and f = r/R, where R is the matching 
radius.Q We then integrate the expression of G 00 — recall Eqs. ( |18D — in order to obtain H 2 , 
getting 



p2 

Mr) = Y 



^(l- e - 3 )+Sf 2 



(45) 



where we have already imposed the regularity conditions at the origin, Eqs. (22) . The 
matching conditions at the spatial hypersurface yield 

a ( e ~~ 1) t-. 6m . 



whence, 



A^ + 2B = -^ + 2A 1 
e R 6 



6m ( e \ „ 6m 1 . , 

A=—\ , B = Ai r. 46 



Finally, using Eqs. 



(G 22 ) mt (r) = A(f 3 - l)/" -B, 



^2 

where A and B have been given above. At the origin 



Gu(0) = G 22 (0) = A 1 + — y—^j . (47) 



6m / 
~RZ\ 

In Sect. [7.1| we will compare, numerically, our results with those in the model of pOf . 
1 For other choices see e.g. Elf . 
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6.4 Ayon— Beato and Garcfa's models 



In a series of papers, ||23||— 1|26|| , some models of regular, electrically charged BHs with an 
energy- momentum tensor of the form of a SSQV were presented — see also [24]. Their 



importance relied in the fact that the sources that give rise to those spacetimes could be 
linked with non-linear electrodynamics (NED), which besides being a theory by itself, may 
be viewed as a low energy limit of string theory or M-theory. Thus, some plausible models 
of regular BHs — that took into account quantum effects in a clearer way than before — were 
put forward. The features of their models are analogous to the case of Dymnikova's model, 
though with a clear interpretation of the source origin. For the sake of brevity, we will focus 
on the model in . 

The choice there was H(r) = (2m/r)[l — tanh (q 2 /2mr)], for any r > 0. Ours will be: 



H(r) 



A 

r 

2m 
r 



1 — tanh ■ 



B 



< r <R, 
R<r, 



(48) 



where A and B are constants to be determined. The matching conditions imply 



B r 



2B 
' R 



(49) 



Defining x 
in 



A/Aq, y = B / B with A = 2m, B = q 2 /2m, that is the values of the model 



26[], we get 



x 



1 

y 



cosh Xy, 1 + e 



-Xy 



2(1 -X)y, 



(50) 



where A = Bq/R = q 2 /2mR. One has here to solve a trascendental equation in order to 
find the appropriate constants of the interior model. The parameter A is the one controlling 
the set of solutions. In classical electrodynamics, A = 1. We see that there is no solution 
in this case. In the context of General Relativity, A = 1 corresponds to the case where the 
exterior metric becomes flat at a spherical surface. But the choice of H int cannot be zero for 
any positive value of r. Therefore the matching is impossible. The same happens for the 
other models in P3[| -|j24|| . In the following section we will see which type of solutions arise 
for different values of A. 



6.5 Bronnikov's model 



In |29[ a model for static, regular, purely magnetically charged BHs with an energy-momentum 
tensor of the type of SSQVs was proposed. Its interest is two-fold. Again the energy- 
momentum content of the objects was directly connected with NED. Second, it turns out that 
those BHs are the only ones based on a Lagrangian formulation of NED with a Maxwellian 
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behaviour in the weak field limit, regardless of the place the weak limit is taken. The exam- 
ple given there was a GNRSS metric with H{r) = (\q m \ 3 ^ 2 / ar)[l — tanh(ai/|g m |/r)] with m 



—the ADM mass — equal to \q m \ 3 ^ 2 /2a, being q m the magnetic charge. It is then obvious 
that the results of the previous subsection are valid now, just by changing q with q m . The 
difference lies in the fact that now one has magnetic fields and also the theory describing 



NED is different to that of Pfl-[E1 



7 Numerical results 

In order to study the approximate values of R for a given object, one needs to assume a 
particular behaviour of the matter and energy inside the source. As of now, there is no 
agreement at this point. However, following several results, see e.g. |15|-p3|, [[4^]-|]47|], the 



geometry of the core may be described by a dS solution. This has been the assumption 
used in most of the works dealing with regularized BHs. Here we will also include two ex- 



amples with a different behaviour and in Sect. [10] we will draw the main lines of a general 
behaviour. In any case, the aim is to choose those physical models which are as consis- 
tent as possible, the dS model being one of them. In this case, at the core we will have 
^00(0) = — Gn(0) = — (^22(0) = — £733(0) = A 2 = const. Nonetheless, there is no present 
agreement about the scale at which regularization could act. A convenient way to handle 
and integrate this indeterminacy is to set A 2 = 10 3s Ai, being s the free parameter that 
governs the renormalization scale. For instance, if s is around 40, we are then considering 
that regularization takes place at Planck scales, and so on. 

Finally, for the exterior region, in accordance to several recent observations |3|, |J , we will 
assume in what follows that e [10 -10 , 2 x 10~ 8 ] erg- cm -3 . An analysis shows, however, 
that the fundamental contribution comes from the quantum gravitational model describing 
the core, and not from the type of quantum vacuum contribution that is assumed for the 
exterior region or near the surface of the body. 

7.1 Two arbitrary powers 

In this numerical analysis we will consider the uncharged case, because there are no observed 
objects that can be associated with static, charged BHs. If there is charge in the source, 
then the interesting situation involves rotation, which might be eventually connected with 
elementary particles (we refer the reader to [pS]). The relation fl3~3|) is (now q = 0) 



whence 




R = RqvM\\—. — ■ — r, (52) 
w 4(N-2)' v ' 
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171 


s = 30 


s = 40 


s = 50 


M Q 


10" 9 


1Q -19 


10 -2B 


1O 3 M 


10" 8 


1Q -18 


1Q -28 


1O 6 M 


10- 7 


lO" 17 


lO" 27 


1O 9 M 


10" 6 


10 -16 


10 -26 



Table 1: R in cm for various astrophysical and galactic objects and different scales of regularization 
(s = 30 corresponds to a GUT's regularization scale, s = 40 to a Planckian one, etc.). In any case 
R/Lfteg is much bigger than 1 (R/L^ cg ~ 10^ 6+s ^ 2 ^). Therefore, all of them are quite far from 
their corresponding regularization scale. 

where we have put m = Mm & , m & being the Sun mass, and R s = ^/24m /(A2 — Ai). 
The last value only depends on the regularization scale and corresponds to the solution 
for a collapsed object of one solar mass in the case of the "lowest powers" model: R s G 
[3 x 10 21 ~ s , 6 x 10 20 ~ s ]cm. For s = 40 we get R s G [3 x 10~ 19 , 2 x 10" 20 ]cm. Yet we see that 
the object has a quantum size very far from Planckian scales, even if s is bigger. In general 
Rs/Lp\ > 10 13 ! Moreover this result is valid for all N, since for any value of A" we have that 
R G [0.6, l]R s vM. It is obvious that, for any astrophysical object, the final properties are 
very similar. Table [l] comprises different massive objects and regularization scales and their 
associated values of R within this model. 



7.2 Israel and Poisson's model 

We have found that the corresponding model within our family must satisfy 

A = (53) 

am 

In this case, A~ 2 = lim^o Coo = A2, so that 

R 3 = -^7^ = ^— (54) 



6 - a,/A 2 q _ Py/Aztfn 

This model clearly depends on the coefficient (3. For instance, in order to obtain a solution, 
we must have (3 2 < 36/(A 2 Lp 1 ). The natural scale of regularization in this model is the 
Planckian one since from the beginning the coefficient a was related to the Planck length. 
Obviously other regularization scales would simply change Lpi by the corresponding scale. 
Using standard values for A2, that use a Planckian regularization scale, and the fact that 
(3 2 should be at most of order unity M, |44j], we get R ~ \f~M x 10~ 20 cm. This result is in 



complete agreement with the foregoing values, even though the models possess very different 
functions H(r). 
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7.3 Dymnikova's model 

From Sect. IO and the assumption of a dS core, we have 



whence 



A2-A1 Ve-2, 

Comparing this result with the one in Eq. ([52]), we get that R = A-Rxwopowers, with .84 < A < 
1.34, for any N. Therefore, R is again of the same order of magnitude, despite the differences 
in the choice of the profile of the energy density and of the tangential pressures. 

Comparing now the model proposed here with the original one in p0| , we see that both 
yield similar conclusions (in the instances they can be compared). For example, in the 
mentioned work, a characteristic radius was found for the collapsed body. Its expression is 



R c = y6m/(A 2 - Ai), what yields me - 2)/(e - 1)R ~ 1.34i?. Besides having a different 
description for this in our model, the values of the coefficients A and B are also quite different 
(numerically). 

7.4 Ayon— Beato and Garcfa's model, and Bronnikov's model 

In the papers dealing with those models, there is no analysis of the orders of magnitude of an 
eventual characteristic radius. The only such condition on these model is to have an event 
horizon. We can now compute which are the ranges of R corresponding to different cases of 
A. 

First of all, Eqs. (|50"D only have solution for < A < 1. Therefore, extremely charged 
objects (those with \q\/m >> 1), cannot be described within the present framework. This 
would require R » m, so that the regularized object would not be a BH but a "visible" ob- 
ject, such as an electron (its size, though, being bigger than the classical radius, or Compton 
size, q 2 /2m). 

For strongly charged objects, i.e. \q\/m ~ 1, we get that, in order to have a BH, 
1/2 < R/m < 2. Thus, the regularized object is of a similar size as that of the event 
horizon. Much bigger than in the uncharged case. 

The solution given by f26|| , i.e. A = 2m, B = q 2 /2m, can only be valid now for very 



weakly charged objects, \q\/m << 1, and satisfying R/m < q 2 /m 2 . They showed that their 
model was acceptable for \q\/m < 1.05. Now, we see that the values of A, B in our model 
change for most of these cases. 

The same is valid for Bronnikov's model, just changing the electric field by a magnetic 
one. Nevertheless, rotation should be introduced in such case — at least when \q\/m is not 
very small. 
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8 Horizons and an interpretation of the regularized 
BH 



Looking at Eq. (30) in and comparing it with our result 

g 00 = -l + H 2 , (57) 

we realize that, substituting here our corresponding H 2 for that model, these expressions 
turn out to be very similar, except for a possible overall sign difference, due to the different 
signatures (e.g. (+,-,-,-) instead of our (-,+,+,+)). We conclude that the same structure 



for the horizons and Cauchy hypersurfaces is obtained. In ||21|| , the solutions are obtained 
by approximations of the exact solution, so that these results and ours are really coincident 
(the relative error with respect to both exact solutions being completely negligible). 

In general, the horizons result from the cancellation of goo- Thus we are left with a 
general set of horizons. A global study for all the candidates encountered, has not yet been 
carried out. We could focus on examples, and try then to extract some general features 
from them, but we do not find this of primary importance. The main point is here, in fact, 
that the matching occurs at a radius which is substantially smaller than the Schwarzschild 
radius of the object. Therefore we will always have a typical exterior, a vacuum transition 
region extending until the matching with the object happens, and a quantum-dominated 
interior, which finally converges to a de Sitter core. In the vacuum interior region and in 
some part of the quantum object, the role of t and r are not interpreted as usual (d t changes 
its character). This is the reason for adequately treating the horizons: to see where exactly 
such changes appear. But, we can still perfectly agree in ordinary physical terms without 
requiring a general determination of the precise radii at which horizons occur. 



Moreover, in ||4 1 , the authors studied the stability of the model. The same considerations 
there hold for our whole family of solutions, as can be easily seen after a careful analysis. 

Finally, there is still the issue of the topology of the solutions, which is connected with 
the possibility of a "universe reborn" in the extended spacetime. Its general structure can 
be found in |4£| for the case where the sources satisfy weak energy conditions (see next 



section). There, it was shown that the topology of any regular BH, satisfying the weak energy 
conditions, should be similar to that of a singularity-free Reissner-Nordstrom spacetime. 
However, there are relevant solutions in our family that violate the weak energy conditions 
(WEC). It would be worth studying what happens in those cases. 



2 With respect to the other models encountered here, we have found that the results are rather similar to 
those in Dymnikova and Soltysek's model pifl. 
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9 Energy conditions 



A common point when dealing with the avoidance of singularities is to show that the energy 
conditions required in the singularity theorems (see e.g. J7|) fail to be valid. 

Here we will study the strong energy conditions (SEC), the weak energy conditions 
(WEC), the null energy conditions (NEC) and the dominant energy conditions (DEC), within 



the GNRSS family (see [35[ for the case of a general spherically symmetric spacetime). The 
SEC are related with the formation of singularities in the collapse of an object. The WEC 
are directly related with the energy density measured by an observer. The NEC are useful 
in order to include some spacetimes which violate the first two, but are predicted by some 
quantum models, e.g. AdS. Finally, the DEC are in fact related with the causal structure of 



the energy-matter content of a spacetime [49].(j Even though an analysis of energy conditions 



helps to understand the physics of a model, one has to be cautious on ascribing to them 
more relevance than they actually have. In several systems, mainly when quantum effects 
play a fundamental role, they all may be violated with less difficulty (see e.g. the review in 

0). 

Let {e a }, a = 0,1,2,3, be a dual vector basis of the cobasis used in |3]^, defined by 
® b e a = 5 b , b = 0, 1, 2, 3. Any timelike vector field, V, in the manifold can be represented by 

V = A b e b , (A°) 2 = l + E(^) 2 , (58) 

i=i 

where A b are some functions. 

On the other hand, from the results of Sect. ||, the Ricci tensor is 

Ricci = .M© ® ©° - 1 ® 1 ) + # 22 (© 2 <8> © 2 + © 3 ® @ 3 ), (59) 
where ® denotes tensor product. A similar expression holds for the Einstein tensor. 



9.1 Strong energy conditions 

SEC require Rvv = R a bV a V b > 0, for all V. From the expressions above, we obtain 
Rvv = Roo + (Roo + -R22)[(^4 2 ) 2 + (^4 3 ) 2 ]- Taking into account the expressions of the Ricci 
and Einstein's tensor given in Sect. |J we get R vv = G 22 + (G 00 + G 2 2)[(^ 2 ) 2 + (^ 3 ) 2 ]- 
Finally, using Einstein's equations, and the fact that A 2 , A 3 are free, we get 

SEC^p + p 2 >0, p 2 >0 (60) 

where p is the energy density measured by e*o, 8ivp = Goo and p 2 is the tangential pressure 
(or stress) of the source, 8iip 2 = G 22 . This is the usual representation of SEC. However, 

3 Let us notice by passing that their Eqs. (2.25), expressing the DEC, are wrong. For our case, the correct 
ones are given in Eqs. (p4|). 
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the GNRSS family allows for a new, and more useful, expression. Indeed, as mentioned 
elsewhere, p 2 = — (p + rp'/2), where ()' = d()/dr. Therefore, we can write 



SEC^j9 2 >0, p'<0. (61) 

In all the examples given before SEC are violated. This is natural since they are regular. 
Particularly, SEC are violated for r < Rsec, with 

Rsec = N Rsec = (-^ 4(6fl"^ ^)) R> Rsec = -GSR, Rsec = -835 

where, all the quantities have been defined in Sect. and the solutions correspond to the 



two-power model, the Israel-Poison's model, Dymnikova's model, and Ayon-Garcia's |26 
and Bronnikov's model, respectively. Indeed, SEC are violated in a main portion of the 
object, i.e. Rsec ~ R- For the evaluation of the Israel-Poison's model, we have used the 



same numerical values as in Sect. |7J]. In the case of Dymnikova's model the value displayed 
corresponds to the case A 2 >> Ai. For any other case with A 2 > A x > or A 2 < < A 1; 
as expected, SEC are violated from bigger values of Rsec- Finally, in the latter case, one 



should evaluate B for different possibilities (see Sect. |6]4] and the next case 



9.2 Weak energy conditions 

Following analogous steps, one finds, for WEC (Gyv > 0, for all V) 

WEC «-> p > 0, pf < 0. (62) 



It turns out that WEC are satisfied in the models of Sects. pTT] , Sect. |0| and Sect. |6T3| , 
very easily for any value of r (e.g. for de Sitter core, p' = 0). One only needs to impose 
Ai < A 2 . 



Let us now consider the series of models in Sect. |6.4| , |6.5| . We have already seen that 
< A < 1. This implies that y, x > 0, and hence that A, B > 0. In general, we have 
H(r) = (A/r)[l-tanh(B/r)). We then get 8vrp = (l/r 2 )(Fr)' = (AB/r 4 ) cosh" 2 (S/r) and 
87rp' = (2AB/r 5 ) cosh -2 (B / r) x [—2 + (B/r) tanh(5/r)]. The energy density p is positive 
for any r, although p' may become positive. To see this, we first solve p' = 0. Its solution is 
r ~ .485. Therefore, we have: For r < .485 fl r < R (outside the body WEC are satisfied), 
WEC are violated. 

In the model of Refs. [26|, ||2~9fl , one has B = q 2 /2m and |g| < 1.05m. This gives that 



WEC are violated for r < .27m, already far away from the core. 

In our revisited model, we have basically two different possibilities. First, for weakly 
charged sources, i.e., those with \q\/m << 1, the conclusions are the same as for the model 
in [2~B\, [29|. Second, for sources with \q\/m ~ 1, we have — recall Sect. |7.4j — m/2 < R < m, 



for a BH. Two limiting alternatives appear. 
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The first one is that R — > 2m. In this case y ~ 1 and, therefore, WEC are violated for 
r < .27m. The other one is that R — > (m/2) + e, with e << 1. Now, y ~ l/4e. WEC are 
violated for r < .06m/enr < 2m, that is everywhere inside the source. In conclusion, WEC 
are again violated almost everywhere. 

Finally, if one lets R > 2m (one does not have now a BH, but a "visible" object) big 
enough to have A < 1 for any |g|, m, we get that WEC are violated everywhere in the object. 

This adds a new (elementary) example to the violation of WEC when quantum effects 
play an important role (see |50| for a recent review) and shows clearly that, although energy 
conditions do help understanding the models, they should not necessarily restrict the search 
for new solutions (Fig. 1). 



9.3 Null energy conditions 

In the case of NEC, V is a null vector field, V ■ V = 0, and requires the evaluation of 
R ab V a V b = G ab V a V b > 0,VV. One obtains 

NEC <-> p' < 0. (63) 

Thus one sees, that a necessary condition common to SEC, WEC and NEC is that the 
energy profile of the sources be a non- increasing function. NEC are satisfied in the models 
of Sects. |BTT] , |£T2] and |6]^, for A 2 > Ai, regardless of the signs of A 2 or Ai. In the models 



of Sects. |6.4| , |6.5| , NEC are violated in the same regions as WEC are, contrary to the belief 
expressed in PD|. 



9.4 Dominant energy conditions 

DEC are satisfied if and only if |T° | > |Tj|, i, j = 1, 2, 3. For the GNRSS family one gets 

DEC <-> signp = sign(— p) = sign(p' + 4p/r). (64) 

Two immediate consequences are, that if p changes its sign, DEC are violated, and if WEC 
are violated in a region with p > 0, then DEC are also violated. 

Let us turn now to the models considered here. In the case of the two-powers model we 
will assume Ai > 0. In this case, p is positive everywhere. WEC were satisfied in these 
models. However, DEC may be violated. A study of the sign of p' + 4p/r tells us that SEC 
are satisfied for f N ~ 2 < 4A 2 / (A 2 — Ai) x 1/(N + 2). Now the question is whether f is less 
than 1 or not. 

Obviously, for any N exceeding N* = 4A 2 /(A 2 - Ai) - 2 = 2(A 2 + Ai)/(A 2 - Ai), we 
have that DEC are violated. One may ask whether this is too odd or easy. Since one expects 
A 2 >> Ai, we readily get A^* ~ 2. This, together with the fact that A^ must be bigger than 2 
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p(r) 




r 



Fig. 1 

Figure 1: Plot of the density p, in arbitrary units, in terms of the coordinate r. RN means 
Reissner-Nordstrom, dS means de Sitter, Schw means Schwarzschild and AdS means anti de Sitter. 
Prn = e 2 /r 4 , p dS = A 2 , pschw = and p A dS = In the region r £ [0,R], p can be any 

(smooth) function matching continuously with p at the center and at the surface. Regions where 
p is increasing violate SEC, WEC and NEC. These are clearly most but not all possibilities. From 
the plot, e 2 /R 4 < 8vrA 2 , if one wants that the model fulfills WEC or NEC. The addition of an 
external A simply shifts the horizontal axis a quantity Ai. 
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(in order to be singularity- free), implies that, in practice, DEC are violated in these models 
— recall WEC are satisfied throughout. 

For the model of Sect. |5^2] , we can assume B, C to be positive (for, if BC < one gets a 
negative Schwarzschild's mass outside the body and if B, C < 0, B, C can be substituted by 
\B\, \C\). Following a similar analysis as with the previous models, one gets that DEC are 
satisfied for r < r*, where r* = (2B /C) 1 / 3 '. Therefore if r* is less than the matching R, there 
is a region where DEC are violated. This is the case of our corrected model. Incidentally, 
in the original model, DEC are violated for r > r*. This conclusion is against physics, since 
far enough one expects Schwarzachild's solution to be valid, and it is a vacuum's solution 
with no problems in its causal structure. Therefore, the corrected version not only describes 
a more realistic picture but also solves this undesired property. 

Turning back to our corrected model we have to answer whether r* may be smaller than 
R — see the expressions given in Sect. IfO. We get that for ft > ft + or ft < ft-, where 



ft + = 6( v / 3 - 1) x ^(L^Aa) and /3_ = -6( v / 3 + 1) x ^(L^A 2 ). On the other hand, it is 
expected, |51[, that LpjA 2 ~ (9(1) and, consequently, ft + , ft- are of order unity. Therefore, 
even though there are several parameters for which DEC may hold, there are also many others 
for which DEC will fail. A more definite answer can only be provided after a particular field 
model is chosen, what will yield a particular value of ft. What is this plausible field model 
remains, as of now, unknown. 

For the next model (the one in Sect. |6.3| ), it is easy to show that DEC are satisfied 
throughout the source if A 2 > Ai > 0, as expected. This is contrary to the other models, 
since this one departs from them through the causal connection in its stress-energy content. It 
is to be noticed that DEC give a new input to understand the models. (The case A 2 < A 1 < 
also satisfies DEC, whereas the rest of possibilities violates them). 

Finally, for the models in Sects. |6]4], |B~5] , as p is positive and p' is positive near the core, 
DEC are violated together with WEC. 

Some concluding remarks are in order. First, although DEC are known to be different 
from WEC, here we see more: it turns out that in cases with p > 0, DEC are more restrictive 
than WEC. Another consequence is that DEC violation and the spacetime region where it 
occurs are not related. That is, DEC may be violated in regions where H(r) is larger or 
smaller than 1. It happens, however, that after substituting expected numerical values for 
the physical parameters involved, the values of H(r) where DEC is violated belong mainly in 
the region where H(r) > 1 and a "signature change in spacetime has occured". The region 
with H(r) < 1 is then at Planckian (regularization) scales and can thus be forgotten. On 
the other hand, when WEC are violated, one usually accepts that the energy-matter content 
of the model can no longer be described by a classical matter source model. However, DEC 
deserve some especial attention. 

These remarks impel us to further interprete the violation of DEC from the causal inter- 
pretation of DEC [[49] . A possibility is that the breakdown of causality in matter interaction 
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may be interpreted in similar — though properly adapted — terms as is the Einstein- Podolski- 
Rosen paradox interpreted in Quantum Mechanics. 

If this is so, or something similar can be proven, DEC may be a more natural sign of 
quantum effects in matter than WEC, for the case of positive densities. This point deserves 
further investigation. 

We will now analyze the main features arising when one replaces the de Sitter core by a 
different spherically symmetric solution. 

10 The matching of static spherically symmetric space- 
times 

In previous sections we have worked with the assumption that the energy-momentum tensor 
satisfies Tq = T±. Now we would like to make the first steps towards the general case where 
Tq and Tl may be independent of each other. Therefore, our aim here is to match two space- 
times that share the existence of an integrable Killing field and spherical symmetry. In order 
to get the most natural junction, we need to take profit of both symmetries exhaustively. 
The metric can always be written, for any of them, as 

ds 2 = g AB (R) dx A dx B + G 2 {R)dtt 2 , (65) 

where A, B = T, R. 8t has been chosen to be the integrable Killing vector and dQ 2 = 
<i# 2 + sin 2 9 dip 2 . Moreover, if G'(R) = 0, we already know that they belong to the generalized 
Nariai family, in which case they only match with another member of its own family. Thus, 
we will only deal with the situation G'(R) ^ 0. In this case, a direct redefinition of the R 
coordinate allows us to write 

ds 2 = g AB (r) dx A dx B + r 2 dtt 2 , (66) 

where A, B = T, r. 

Spherical symmetry has thus been completely used. We now extract consequences from 
the presence of dr- The natural thing to do is to identify both vector fields, i.e. &n = 9t 2 - 
However this is not a right choice, in general, because if a Killing vector is multiplied by 
a constant factor, the resulting vector field is obviously a Killing vector field. Therefore, 
normalizing each Killing vector, when possible, gives the natural way to identify them. This 
is indeed implemented in the junction process, if the hypersurfaces are spacelike or timelike 
everywhere. On the contrary, in the general case, we cannot rely on such normalization. 

Any metric of interest (to our purposes) can be written as (recall the coordinate change 
to obtain @, setting now F = 1 — H) 

ds 2 = -(1 - H) dT 2 + — dT dr + ] -^JL dr i + r 2 ( 67 ) 
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where H, g ^ are functions of r only. Looking back to expressions (|12"D and the coordinate 
changes mentioned there, we will put now dT = g dt + (g — gQ^dr, where g is a constant, 
that will be related with the function g, as we shall see in a moment. With this coordinate 
change the metric takes the form 

ds 2 = -glil - H) dt 2 + 2G dt dr + F dr 2 + r 2 dQ 2 , (68) 

where G = g 2 (H - 1) + 1 + H(g - g)/g, and F = 2 + g 2 (H - 1) + (g Q - g)[2H/g + 2/g 2 g + 
(g -g)(H+l)/g 2 g). 

The junction conditions (for any type of hypersurface, see ||52|| ) are then 



[r] = 0, [t] = (69) 
\H}i 2 + 2[G]if + Ff 2 = 0, (70) 
[F]if - [G] (it - ff) + [H]fi + [G'}f 3 + [H' - (F'/2)}f 2 i - [H']i 3 /2 = 0, (71) 



where [/] = f% — /i, and where we have put H = g^H — 1) + 1. In (J70) and ([71] ) t and 
are either £i, fi or £2, ^2, and A' = dA(r)/dr| r=r m. The same conditions lead, in general, 
to a second order ordinary differential equation for r. In principle there is the possibility for 
asymptotic stopping solutions, i.e. solutions for which r — > const, as t — > 00, and also for 
null ones. The special case r\ = r 2 = R = const, is of great interest, since it constitutes the 
solution towards which any transitory solution should converge. Under this restriction, the 
conditions become, simply, 

[t] = 0, [H] = 0, 2[G]i - [H'}i 2 = 0, (72) 

where t is either t\ or t 2 . Choosing g as g% one gets [G] = (the same result comes out 
directly in the case when the normal vector of S is non-null). The last conditions become 
then [H'] = 0. Thus, the conditions emerging from the matching of two spherically symmetric 
spacetimes with an integrable Killing vector field are, for the case r = R— const, and taking 
the maximum identification between them, 

[H] = 0, [H'\ = 0, (73) 

where H = g^(H — 1) + 1. An intrinsic characterization, valid for any representation of the 
form (|65|) or (|66|) (the ones most often dealt with in the literature) is H = — • ^) + 1, 
g s = [G'/\ det(gAB)| 1,/2 ]r=i?, where ^ is the Killing vector associated with the staticity of the 



r 



solution (in some regions) of (65) or (66). Finally, notice that the first condition on H is 
nothing but the requirement of the mass function to be continuous across the hypersurface, 
while the second one is related with the continuity of the radial stress, or pressure (see 
e.g. (0)). Needless to say, if one restricts oneself to the family of metrics in fll2"D , one gets 



the conditions of Sect. 4.2 



25 



11 An application to supersymmetric stringy black 
holes 



The semiclassical expressions for supersymmetric stringy black holes are well-established (see 
e.g. Q and references therein). There are also other objects of interest, such as black 
strings, higher dimensional black holes, etc. In all cases, one looks for a correspondence 
principle with general relativistic black holes. This transition is usually reflected in the 



strength of the coupling constant, or the entropy (see e.g. |53]-|56] and references therein). 
Here we take a complementary viewpoint. 



The most interesting case to our aims is that of a self-gravitating string (see e.g. [5^, |56|j). 
However the necessary ingredients — specially the corresponding spacetime metric — in order 
to tackle this problem are still under study. Here we will consider the most simple (and widely 
considered) case, that of a supersymmetric back hole. 

A family of such black holes, related with electrically charged black holes, is given by 
(see |53|, |54| for details) 



ds 2 



r 



el/2. 



(l-^) 'dr 2 + r 2 dn 2 



(74) 



where f(r) = Ilf=i[l + ( r o sinh 2 atj/r)], and where the ai are related with the integer charges 
of the D-branes being used. If the correspondence occurs at a constant value of r, we get 



r\ + ro sinh a 
2m _ Q 2 



r 2 / 2 V4 (r 2 ; 



^-1 



= R — const. 



2m 2Q 2 
~~R 2+ ~R T 



1 + 



If) L2f\ r) r 2 \) 



2/ 



(75) 
(76) 

(77) 



where we have used g% = G f \ r =G-i(R), G(r) = r/ 1 / 4 (r), and ^ • £ = — /~ 1 ^ 2 (1 — r/r ). The 
subscript S 2 means that all these quantities refer to the interior region, to be evaluated at 
r = r 2 . For the exterior metric, we have put Q!j = atj = a, for all i, j, because the exterior 
metric is that of a Reissner-Nordstrom black hole, for which 



2m = r cosh 2a, Q 2 = Tq sinh 2 a cosh 2 a, 

ro = 2\/m 1 — Q 2 , 2 sinh 2 a = — 1 + m/y/m 1 — Q' 2 , 



(78) 



where m is the (ADM) mass of the black hole and Q is its electric charge. Since /2( r 2 ) = 
rd=i[l + ( r o sm h 2 o;j/r)] S2 , the above conditions yield R as a function of six of the seven 
parameters, M, Q, (r )2, on- Detailed analysis shows that these conditions are easily fulfilled 
when r — > 0, «j — > ±oo, with r sinh 2 aj fixed. The resulting R is very close to -Ro = 
m + y/m 2 — Q 2 , i.e. the event horizon of the black hole. We remark that r/ 1//4 (r) is the 
radial coordinate which has a direct interpretation in terms of the "size" of the object, and 
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not r alone. All this being in complete agreement with the expected transitions for extreme, 
and nearly extreme, supersymmetric black holes. The same idea should be extended to self- 
gravitating strings when their (4-dimensional) spacetime metric is obtained. For instance, 
the expected order of magnitude of R found in |55[] , should be recovered. This issue will be 
the matter of subsequent research. 



12 Final remarks 

The first thing to be noticed is the intrinsic freedom present in our model, which is as large as 
the measure of the set of analytic functions of one variable. This is a very rewarding feature, 
since it allows to impose further restrictions coming from new, more accurate proposals. In 
particular, it will be a helpful tool when trying to find explicitly a quantum field responsible 
for the Gn and G22 in the fundamental uncharged case. For comparison, in all previous 
works, based on individual models, the prospective of finding a quantum field related with 
their energy-matter content was hopeless. To that end, we would like to draw the attention 



to | 57fl , where a useful framework to deal with the interior region is given. 



In the charged case, let us notice that any GNRSS spacetime can be linked with a solution 



to NED (see [5S]). Of course, the case of Schwarzschild solution is a solution with zero charge 
and Reissner- Nordstrom one, the only one which is linear, i.e., Maxwellian. Therefore, the 
whole family of GNRSS metrics has indeed an immediate interpretation in terms of a field 
theory which is well established when the object is electrically or magnetically charged. This 
is another useful result. A carefull study of this fact will be reported elsewhere. 

Finally, one can free the requirement that there must be an event horizon. The objects 
would then become "visible" and the study of the entropy of the solutions as well as their 
associated Hawking temperature would bring some clues on the time evolution of (classically) 
static black holes. 



13 Conclusions 

In this work we have investigated, under quite general conditions, the question of using Ein- 
stein's theory of gravitation — extended to include semiclassical effects — with the purpose to 
constraint the physical structure of the emerging spacetime solutions that might be suitable 
for the description of the interiors of non-rotating black holes.0 

In the first part of the work we have exploited the idea that vacuum polarization may 
indeed play an essential role in the interior region. We have obtained the result that only two 

4 The rotating case, which is of major astrophysical interest, and the rotating and electrically charged one, 
which may be associated with spinning particles, seem to yield results very similar to the ones presented 
here, see g§, 
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families fulfill the imposed requirement and, moreover, we have shown that only one of them 
is suitable for representing static black hole interiors, what is certainly a most remarkable 
result. 

Then we have turned our attention to other sources for the core. Given that a promising 
alternative — self-gravitating strings — needs still to be studied in more detail, we have 
started this program by first giving the general conditions to be fulfilled by any spacetime 
with spherical symmetry and having some static region. Finally, we have applied the results 
obtained to a supersymmetric black hole, as a preliminar case. We have seen that, in such 
situation, the matching is generically compatible, including the case of extreme black holes. 
This last setting is precisely the same for which the correspondence between semiclassical 



black holes and stringy ones has been recently confirmed in the literature (see e.g. [55, pq|). 

Briefly, our overall conclusion is the following. First, the results in the first sections have 
opened a new window for the search of a compatible quantum field that, once regularized, 
may yield the same result for, at least, a particular energy- momentum tensor inside the 



general family of models considered (for instance within NED, see also ]5"9"|). Second, once 
a corresponding Einsteinian metric associated with a quantum model is known, the scheme 
developed here has been proven to be well suited to check the consistency of the involved 
physical parameters and even, in some cases, to assign explicit values to them. 
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